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Abstract. In this paper we present a new methodology for solving multiobjective integer linear programs using 
tools from algebraic geometry. We introduce the concept of partial Grobncr basis for a family of multiobjective 
programs where the right-hand side varies. This new structure extends the notion of Grobncr basis for the 
single objective case, to the case of multiple objectives, i.e., a partial ordering instead of a total ordering over 
the feasible vectors. The main property of these bases is that the partial reduction of the integer elements in 
the kernel of the constraint matrix by the different blocks of the basis is zero. It allows us to prove that this 
now construction is a test family for a family of multiobjective programs. An algorithm 'a la Buchbergor' is 
developed to compute partial Grobncr bases and two different approaches are derived, using this methodology, 
for computing the entire set of efficient solutions of any multiobjective integer linear problem (MOILP). Some 
examples illustrate the application of the algorithms and computational experiments are reported on several 
families of problems. 



1. Introduction 

The multiobjective paradigm appeared in economic theory in the nineteenth century in the seminal works 
by Edgeworth [Mj and Pareto [34j to define an economic equilibrium. Mathematically, the multiobjective 
optimization approach consists of determining the maximal (minimal) elements of a partially ordered set. This 
problem was already addressed by Cantor [7] , Cayley ^ and Hausdorff [25] at the end of the nineteenth century. 
Since then, multiobjective programming (including multicriteria optimization) has been a fruitful research field 
within the areas of applied mathematics, operations research, and economic theory. Excellent textbooks and 
survey papers are available in the literature, the interested reader is referred to the books by Sawaragi, Nakayama 
and Tanino [Mj, Chankong and Haimes [H], Yu [SU], Miettinen [3S] or Ehrgott, Figueira and Gandiblcux [20], 
and to the surveys in [T7] and [T^ . 

The importance of multiobjective optimization is not only due to its theoretical implications but also to its 
many applications. Witnesses of that are the large number of real-world decision problems that appear in the 
literature formulated as multiobjective programs. Examples of them are flowshop scheduling (see |29|). analysis 
in finance (see jl7| . Chapter 20), railway network infrastructure capacity (see |13j). vehicle routing problems 
(see [301 [5H]) or trajectory optimization (see [H]) among many others. 

Multiobjective programs are formulated as optimization (without lost of generality, we restrict ourselves to 
the minimization case) problems over feasible regions with at least two objective functions. Usually, it is not 
possible to minimize all the objective functions simultaneously since the objective functions induce a partial 
order over the vectors in the feasible region, so a different notion of solution is needed. A feasible vector is 
said to be Pareto-optimal (efficient or non-dominated) if no other feasible vector has componentwise smaller 
objective values, with at least one strict inequality. 
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This paper studies multiobjective integer linear programs (MOILP). Thus, we assume that aU objective functions 
and constraints that define the feasible region are linear, and that the feasible vectors have non-negative integer 
components. 

There are nowadays several exact methods to solve MOILP (see [I7|). Two of them claimed to be of general 
use and have attracted the attention of researchers over the years: multiobjective implicit enumeration (see 
[51j . [52j ) and multiobjective dynamic programming (see [31j). Nevertheless, although in principle they may be 
applied to any number of objectives, one can mainly find, in the literature, applications to bicriteria problems. 
On the other hand, there are several methods that apply to bicriteria problems but that do not extend to the 
general case. Thus, one can see that there are two thresholds in multiobjective programming, a first step from 
1 to 2 objectives and a second, and deeper one. from 2 to more than two objectives. Thus, most of the times, 
algorithms to solve multiobjective integer problems are designed to compute only the solutions for the bicriteria 
case. Moreover, some methods even do not provide the entire set of Pareto-optimal solutions, but the supported 
ones (those that can be obtained as solutions of linearly scalarized programs). 

It is worth noting that most MOILP problems are NP-hard and intractable (see [16] for further details). Even 
in most cases where the single-objective problem is polynomially solvable the multiobjective version becomes 
NP-hard. This is the case of spanning tree problems and min-cost flow problems, among others (see [51] and 
jl5p. Therefore, computational efficiency is not an issue when analyzing MOILP. The important point is to 
develop tools that can handle these problems and that give insights into their intrinsic nature. The goal of 
this paper is to present a new general methodology for solving MOILP using tools borrowed from algebraic 
geometry. The usage of algebraic geometry tools in integer programming (single criterion) is not new (see [lOj , 
|26j . [46] . [27] . [49] . [48]). The main idea is to compute a Grobner basis for certain toric ideals (related to the 
constraints matrix) with a monomial order induced by the objective function. 

Grobner bases were introduced by Bruno Buchberger in 1965 in his PhD Thesis [6|. He named them Grobner 
bases paying tribute to his advisor Wolfgang Grobner. This theory emerged as a generalization, from the one 
variable case to the multivariate polynomial case, of the greatest common divisor in an ideal sense. One of 
the outcomes of Grobner bases theory was its application to Integer Programming, firstly published by Conti 
and Traverso [10] . After this paper, a number of publications using Grobner bases to solve integers programs 
appeared in the literature. 

In [26| , Hosten and Sturmfcls gave two ways to implement Conti and Traverso algorithm that improve in many 
cases branch-and-bound algorithm to solve, exactly, integer programs. Thomas presented in [IS] a geometric 
point of view of the Buchberger algorithm as a method to obtain solutions of an integer program. Later, Thomas 
and Weissmantel [48] improved the Buchberger algorithm in its application to solve integer programs introducing 
truncated Grobner bases. At the same time, Urbaniak et al [H] published a clear geometric interpretation of 
the reduction steps of this kind of algorithms in the original space (decision space). The interested reader can 
find excellent descriptions of this methodology in the books by Adams and Loustanau [5] , Sturmfels [12] , Cox et 
al [12] or Bertsimas and Weissmantel [5], and in the papers by Aardal et al. [T], Sturmfels [43], [44], Sturmfels 
and Thomas [15] and Thomas 07] 

The main contribution of this paper is to adapt some tools from algebraic geometry to solve multiobjective 
integer linear programs. We present in this paper an algorithm to solve exactly multiobjective problems, 
i.e. providing the whole set of Pareto-optimal solutions (supported and non-supported ones). One of the main 
advantages of our approach is that the number of objective functions does not increase significantly the difhculty. 
A new geometric approach of the concept of reduction based on a partial ordering is given. This reduction allows 
us to extend the concept of Grobner basis when a partial ordering rather than a total order is considered over 
N". We call these new structures partial Grobner bases or p-Grobner bases. We prove that p-Grobner bases 
can be generated by a variation of the Buchberger algorithm in a finite number of steps. The main property of 
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a p-Grobncr basis is that, for each pair in Z" x Z!f: with first component in Ker{A), the reduction by maximal 
chains in the basis is the zero set. 

Wc propose two different algorithms to solve multiobjcctive integer programs based on this new construction. 
Our first algorithm consists of three stages. The first one only uses the constraint matrix of the problem and it 
produces a system of generators for the toric ideal I a (or its geometric representation, 3a)- In the second step, 
a p-Grobner basis is built using the initial basis given by the system of generators computed in the first step. 
This step requires to fix the objective matrix since it induces the partial order used in the reduction steps. Once 
the right-hand side is fixed, in the third step the Pareto-optimal solutions are obtained. This computation uses 
the new concept of partial reduction of an initial feasible solution by the p-Grobner basis. 

This algorithm extends, to some extent, Hosten-Sturmfels' algorithm |26j for integer programs because if we ap- 
ply our method to single-objective problems, partial reductions and p-Grobner bases coincide with the standard 
notions of reductions and Grobner bases, respectively. 

Our second algorithm is based on the original idea by Conti and Traverse [T^. It consists of using the big- 
M method that results in an increasing number of variables, in order to have an initial system of generators. 
Moreover, this approach also provides an initial feasible solution. Therefore, the first step in the above algorithm 
can be ignored and the third step is highly simplified. In any case, our first algorithm (the one extending Hosten- 
Sturmfels approach) has proved to be more efficient than this second one since computation of a p-Grobner 
basis is highly sensitive to the number of variables. 

Both algorithms have been implemented in MAPLE 10. In this paper we report on some computational 
experiments based on two different families of problems with different number of objective functions. 

This paper is organized as follows. In Section [2] we give the notation, the formulation of the problem, and its 
algebraic codification. In this section we also introduce the notion of test family and its geometric description. 
Section [3] presents the definition of p-Grobner basis, based on the notion of partial reduction. Here, we also 
state the relationship between test families and p-Grobner bases: the reduced p-Grobner basis for a family 
of multiobjcctive programs varying the right-hand side coincides with the minimal test family for that family. 
At the end of the section, an illustrative example is presented. Section |4] is devoted to present the results 
of the computational experiments and its analysis. Here, we solve several families of MOILP, report on the 
performance of the algorithms and draw some conclusions on their results and their implications. 



2. The problem and its translation 

The goal of this paper is to solve the multiobjcctive integer linear program (MOILP) in its standard form: 

min {c^ X, . . . , x) 
s.t. 

n 

(1) ^ Xj =h i= 1,. ..,m 

Xj eZ+ j = l,...,n 

with bi nonnegative integers, Xi non negative and the constraints are defining a polytope (bounded). Let us 
denote by ^ = (ay) e Z™''", b = {bi) G Zip and C = (c^) e In the following. Problem ^ will be 

referred to as MIPa,c{^) and we denote by MIPa.c the family of multiobjcctive problems where the right-hand 
side varies. 



4 



VfCTOR BLANCO AND JUSTO PUERTO 



The reader may note that there is no loss of generahty in our approach to multiobjective integer hnear program- 
ming since any general multiobjective integer linear problem with inequality constraints and rational components 
in A, b and C can be transformed to a problem in the above standard form. 

It is clear that the problem MIPA,c{b) is not an usual optimization problem since the objective function is a 
vector, thus inducing a partial order among its feasible solutions. Hence, solving the above problem requires an 
alternative concept of solution, namely the set of non-dominated or Pareto-optimal points (vectors). 

A feasible vector x e R" is said to be a Pareto-optimal solution of MIPa,c{^) if there is no other feasible vector 
y such that 

Cjy < CjX = 1 , . . . , fc 

with at least one strict inequality for some j. 

If x is a Pareto-optimal solution, the vector (cix, . . . ,Ck x) is called efficient. 

We say that a feasible point, y, is dominated by a feasible point x if Ci x < Ci y for alH = 1, . . . , fc, with at least 
one strict inequality. According to the above concept, solving a multiobjective problem consists of finding its 
entire set of Pareto-optimal solutions, including those that have the same objective values. 

From the objective function C, we obtain a partial order over Z" as follows: 



x<c y ■< 



X -<c y :4==> C X y or x = y 
where Cx ^ Cy stands for Cx < Cy and Cx ^ Cy. 

Notice that since C G Z™^", the above relation is not complete. Hence, there may exist incomparable vectors 
(those X, y G Z" such that neither x -<c y or y <c x)- We use this partial order, induced by the objective 
function of Problem MIPa,c a-s the input for the multiobjective integer programming algorithm developed in 
this paper. 

Remark 2.1. The above order distinguishes solutions with the same objective values and handles them as in- 
comparable. This order can be refined so that those solutions with the same objective values are not incomparable. 
Consider the binary relation: 

C X ^ C y or 
Cx = Cy and x -<iex y 

This alternative order allows us to rank those solutions that have the same objective values using the lexico- 
graphical order of their components. 

Let us consider the following equivalence relation in Z"; 

X V :<S=^ Cx = Cy 

The above partial order, :<c, allows us to solve a simplified version of the multiobjective problem. In this 
version, we obtain solutions in Z"/ where x y :<^=^ Cx = Cy. The reader may note that when solving 
the problem with the order <c, one would obtain only a representative element of each class of Pareto-optimal 
solutions (the lexicographically smallest). With those efficient values, {wi, . . . ,ft}, the remaining solutions can 
be obtained solving the following system of diophantine equations, in x, for each Vi, i ~ 1, . . . ,t: 

Cx — Vi 
Ax =b 

X e Z 'l 

Remark 2.2. In some cases, the order -<c can be refined to be adapted to specific problems. This is the case 
when slack variables appear in mathematical programs. Two feasible solutions {x,si) and {x,S2), where si 
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and S2 are the slack components, have the same objective values. The order -<c considers both solutions as 
incomparable, although they are the same because we are looking just for the x-part of the solution. In these 
cases, we consider the following refined partial order in Z" x Z'', 

I Cx = Cy and s ^lex s 

where x,y E Z!f_ are the actual decision variables and s,s E Z!j_ the slack variables of our problem. 



In the following wc will use the partial order unless it is explicitly specified. 

Our matrix A is encoded in the set 

Ia {{u,v} ■.u,vE N", u -vE Ker{A)}. 

Let TT : N" — > Z" denote the map x i-^ Ax. Given a right-hand side vector 6 in Z", the set of feasible solutions 
to MIPA,c{b) constitutes 7T~^{b), the preimage of b under this map. In the rest of this paper, we identify the 
discrete set of points 7r~^(6) with its convex hull and we call it the 6-fiber of MIPa.c- Thus, 7r~^(6) or the 
6-fiber of MIPa,c is the polyhedron that is the convex hull of all feasible solutions to MIPA.c{b). 

For any pair {u, w}, with u,v E N", we define the set setlm{u, v) as follows: 

{{u} if V u 
{v} if u <c V 
{u,v} if u and v are incomparable by -<c 

The reader may note that setlm{u, v) is the set of degrees of the leading monomials according to the identification 
{u, v} 1-^ x" — x^ E R[a;i, . . . , Xn], induced by the partial order ^c- 

From the above definition, setlm{u,v) may have more than one leading term, since is only a partial order. 
To account for all this information we denote by J'{u, v) the set of triplets 

J-{u,v) — {(u,v,w) : w E setZm(u, w)}. 

The above concept extends to any finite set of pairs of vectors in N", accordingly. For a pair of sets u = 
. . . , Ut} and v = {vi, . . . , vt} the corresponding set of ordered pairs is: 

J-{u, v) ~ {{ui, Vi,w) : w E setlm{ui, Vi),i = 1, . . . ,t}. 

JF(u, v) can be partially ordered based on the third component of its elements. Therefore, we can see T{u, v) 
as a directed graph G{E, V) where V is identified with the elements of J^{u, v) and ((ui, Vi, w), {uj, Vj, w )) E E 
if {ui,Vi,w), {uj,Vj,w ) E V and w -<c w. We are interested in the maximal ordered chains of G. Note that 
they can be efficiently computed by different methods, see e.g. [4], [37] . 

The above concepts are clarified in the next example. 



Example 2.1. Let u = {(2, 3),(0, 2), (3, 0),(2, !),(!, 1)}, v = {(1, 4),(1, 3), (4, 2),(1, 2),(1, 0)} and <c the partial 
order induced by the matrix 

2 1 

3 5 
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((3,0), (4, 2), (4, 2)) ((2,3), (1,4), (1,4)) 




((0,2), (1,3), (1,3)) 




((1,1),(1,0),(1,1)) 



Figure 1 . Hassc diagram of the graph associated with the data in Example 12.11 

then, setZm((2,3),(l,4)) = {(2, 3),(1, 4)}, setZm((0, 2), (1, 3)) = {(1,3)}, setZm((3, 0), (4, 2)) = {(4,2)}, 
set/m((2, 1), (1,2)) = {(2, 1),(1, 2)} and setlm{{l, I), (1,0)) = {(1,1)}. Now, by definition we have: 

.F(u,v) ={ ((2,3),(1,4),(2,3)),((2,3),(1,4),(1,4)), ((0, 2), (1, 3), (1, 3)), 

((3, 0), (4, 2), (4, 2)), ((2, 1), (1, 2), (2, 1)), ((2, 1), (1, 2), (1, 2)), 

((1,1),(1,0),(1,1))}. 

Figure[^ corresponds to the directed graph associated with .7-'(u, v), according to the partial ordering induced by 
C . There are four maximal chains: 

Mi = {((3,0),(4,2),(4,2)),((2,3),(l,4),(2,3)),((0,2),(l,3),(l,3)),((2,l),(l,2),(2,l)),((l,l),(l,0),(l,l))} 
M2 = {((3,0), (4,2), (4,2)), ((2,3), (1,4), (2,3)), ((0,2), (1,3), (1,3)), ((2,1), (1,2), (1,2)), ((1,1), (1,0), (1,1))} 
Af3 = {((2,3), (1,4), (1,4)), ((0,2), (1,3), (1,3)), ((2,1), (1,2), (2,1)), ((1,1), (1,0), (1,1))} 
M4 = {((2,3),(1,4),(1,4)),((0,2),(1,3),(1,3)),((2,1),(1,2),(1,2)),((1,1),(1,0),(1,1))}. 

For any pair of sets u — {ui, . . . ,ut} and v ~ {vi, . . . ,vt\ with {m;, Vi} G I a, the corresponding set T{u, v) 
may be seen as a set of pairs in Z" x Z" through the fohowing map 

(j) :N" X N" X N" — >Z" x Z\ 
(u,v,w) ^(u — v,w). 

Then, the maximal chains, Fi, . . . , Ft, of the image of ^(u, v) under (j) with respect to the order over the 
second components, clearly satisfy the following properties: 

(1) Fi is totally ordered by the second components of its images via <j> with respect to -<c, for i = 1, . . . , i. 

(2) For all (a,/3) e F„ i 1, . . . ,t, A {(3 - a) = A p. 

The application (p and the above properties allow us to define the notion of test family for MIPa,c- This notion 
is analogous to the concept of test set for a family of single objective integer programs when we have a partial 
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order rather than a total order over N" (see [H]). Test famihes are instrumental for finding the Pareto-optimal 
set of each member M/Pa,c(^) of the family of multiobjeetive integer linear programs. 

Definition 2.1 (Test Family). A finite collection Q = {Qq, . . . ,Gc} of sets in Z" x Z!f. is a test family for 
MIPa^c if only if: 

(1) Q^^ is totally ordered by the second component with respect to -<c, for j = 1, . . . ,r. 

(2) Forallig,h)eg'c,j^l,...,r, Aih-g)^Ah. 

(3) If X G N" is a dominated solution for MIPA,c{b), with h G Z!f., there is some in the collection and 
{g, h) € Q^(j, such that x ~ g -<c x. 

(4) If X G N" is a Pareto-optimal solution for M IPa.c {b), with b € Z" , then for all (g, h) € Q-^j and for all 
j ~ 1, . . . ,n either x — g is infeasible or x — g does not compare with x. 

Given a test family for MIPa.c there is a natural approach to find the entire Pareto-optimal set. Suppose we 
wish to solve MIPA,c{b) for which x* is a feasible solution. 

If X* is dominated then there is some j and (g, h) G Qq such that a;* — 17 is feasible and x* — g -<c x* , whereas 
for the remaining chains there may exist some [g, h) such that x* — g is feasible but incomparable with x* . We 
keep tracks of all of them. 

If X* is non-dominated, we have to keep it as an element in our current solution set. Then, reducing x* by 
the chains in the test family we can only obtain either incomparable feasible solutions, that we maintain in our 
structure, or infeasible solutions that are discarded. 

The above two cases lead us to generate the following set. From x* we compute the set of incumbent solutions: 

IS(x* ) :^ {y* : y* = X* — gj., [gj. , hj. ) is the largest element (17, h) in the chain 
such that x* — g is feasible , i = 1, . . . ,r}. 

Now, the scheme proceeds recursively on each element of the set IS{x*). Finiteness of the above scheme is clear 
since we are generating a search tree with bounded depth (cardinality of the test family) and bounded width, 
each element in the tree has at most r (number of chains) followers. Correctness of this approach is ensured 
since any pair of non-dominated solutions must be connected by a reduction chain through elements in the test 
family (see Theorem 12.11 and Corollarv l2.ip . 

The above approach assumes that a feasible solution to MIPA,c{b) is known (thus implying that the problem 
is feasible). Methods to detect infeasibility and to get an initial feasible solution are connected to solving 
diophantine systems of linear equations, the interested reader is referred to [35j . for further details. 

The following lemmas help us in describing the geometric structure of a test family for multiobjeetive integer 
linear problems. 

Lemma 2.1 (Gordan-Dickson Lemma, Theorem 5 in [TT]). //PC N", P ^ 0, then there exists a minimal 
subset {pi, . . . ,Pm} C P that is finite and unique such that p E P implies pj < p (component- wise) for at least 
one j ~ 1, . . . ,m. 

Lemma 2.2. There exists a unique, minimal, finite set of vectors ai, . . . , ak G N" such that the set Cc of all 
dominated solutions in all fibers of MIPa.c a subset of N" of the form 

k 

i=i 

Proof. The set of dominated solutions of all problems MIPa,c is: 

= {a G N" : 3/3 G N" with A(3 = Aa and (3 a}. 
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Let a be an clement in Cc and /3 a Pareto-optimal point in the fiber n^^^Aa) that satisfies /3 -<c ce- Then, for 
any 7 G N", A{a + 7) = A{(3 + 7), a + 7, /? + 7 e N" and /3 + 7 -<c' + 7j because the cost matrix, C, has only 
nonnegative coefficients. Therefore, a + 7 is a feasible solution dominated by /? + 7 in the fiber n~^{A(a + 7)). 
Then, a + 7 G Cc for all 7 G N", so, a + N" C Cc- By Lemma \2A] we conclude that there exists a minimal set 
of elements ai, . . . , a*; G N" such that Cc = iJj^Mj + N"). □ 

Once the elements ai,...,ak that generates Cc (in the sense of the above result) have been obtained, one 
can compute the maximal chains of the set {ai, . . . ,ak} with respect to the partial order -<c- We denote by 
Cq, . . . ,Cq these maximal chains and set £^ = ljf=i('^t + ^")' where a\ G for i = 1, . . . , fc^ and i = 1, . . . , ^. 
For details about maximal chains, upper bounds for its cardinality and algorithms to compute them for a 
partially ordered set, the reader is refereed to [4]. 

It is clear that, with this construction, we have: Cc — [J C}j. 

i=l 

We now describe a finite family of sets Q-^^ C Ker{A)r\'E^ and prove that it is indeed a test family for MIPa.c- 
Let = {g\c}i^i, where 

(2) ^ ) = ("j - P% ' "j)' = 1, • • • fcj, /c = 1, . . . , my }, 1 = 1, . . . , /i, 

are the maximal chains of Q-^^ (with respect to the order ^c over the second components) and where . . . ,a\, 
are the unique minimal elements of C^^^ and 13}^ , • ■ • , /?™'^ the Pareto-optimal solutions to the problem AZ/Pa^c' ( AaJ ) . 

In the next section we give an algorithm that explicitly constructs Q^^. Notice that for fixed i,j and fc, 

g^j ~ (a*- — /3j^,) is a point in the subspace S ^ {x E Q" : Ax = 0}, i.e., in the 0-fiber of MIPa,c- Geometrically 

— ^ > 

we think of (aj — /3^j,a'j) as the oriented vector g = [aJ■,/3J^,■] in the ^aj-fiber of MIPa,c directed to the 
Pareto-optimal solution /3,^-. The vector is directed from the non-optimal point a* , to the Pareto-optimal 
point /3*j due to the minimization criterion in MIPa,c which requires us to move away from expensive points. 
Subtracting the point = a* — to the feasible solution 7 gives the new solution 7 — -I- /3*j- which is 
equivalent to translating "g'lj by a nonnegative integer vector. 

Consider an arbitrary fiber of MIPa.c and a feasible lattice point 7 in this fiber. For each vector "^f^- in Q^^ , 
check whether 7 — g*^ is in N". At 7 draw all such possible translations of vectors from 5^^- The head of the 
translated vector is also incident at a feasible point in the same fiber as 7 since g^^ is in the 0-fiber of MIPa,c- 
We do this construction for all feasible points in all fibers of AIIPa,c- From Lemma [2.21 and the definition of 
C/^P, it follows that no vector in can be translated by a in N" such that its tail meets a Pareto-optimal 
solution on a fiber unless the obtained vector is incomparable with the Pareto-optimal point. 

Theorem 2.1. The above construction builds a connected directed graph in every fiber of MIPa.c- The nodes 
of the graph are all the lattice points in the fiber and (7,7 ) is an edge of the directed graph if j = 7 — 5,%,- for 
some i, j and k. For each maximal chain in the b-fiber of AH Pa. c, ^ts directed graph has a unique final node 
at each Pareto-optimal solution for MIPA,c{b)- 

Proof. Pick a fiber of MIPa,c and at each feasible lattice point construct all possible translations of the vector 
~g'lj from the set Q]^^ as described above. Let a be a lattice point in this fiber. By Lemma a = a* + for 
some i G {1, • . • , i} and v G . Now, since a^, = (3^^ + v also lies in this fiber, then a^. -<c oi or ct^ and a are 
incomparable. Therefore, ~g\j translated by 1/ G N" is an edge of this graph and we can move along it from a 
to a point ol in the same fiber, such that a -<c a or a and a are incomparable. This proves that from every 
dominated point in the fiber we can reach an improved or incomparable point (with respect to -<c') in the same 
fiber by moving along an edge of the graph. 
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By the construction above, the outdegree of any terminal element in any maximal chain is 0. Therefore, any 
directed maximal path from a dominated point must end exactly at one Pareto-optimal point. □ 

We call the graph in the &-fiber of MIPa,c built from elements in the -<c-skeleton of that fiber. 

The reader may note that from each dominated solution a, one can easily build paths to its comparable Pareto- 
optimal solutions subtracting elements in G-j^- Indeed, let be a minimal element of Cc such that a = 01^ + 7, 
with 7 S N", and let Pi be the Pareto-optimal solution in the Aai-fiher that is comparable with a; and such 
that (3i + J is comparable with (3. Then a = /?i -I- 7 is a solution in the Aa-fiber with /3 -<c a <c oi- Now, 
one repeats this process but starting with a and /3, until a ~ (3. Moreover, the case where a and (3 are 
incomparable reduces to the previous one by finding a path from a to any intermediate point (3 that compares 
with (3. This analysis leads us to the following result. 

Corollary 2.1. In the -<c -skeleton of a fiber there exists a directed path from every feasible point a to each 
Pareto-optimal point, (3, in the same fiber. The vectors of objective function values of successive points in the 
path do not increase componentwise from a to f3. 

Corollary 2.2. The family Q^q is the unique minimal test family for AIIPa^c- H depends only on the matrix 
A and the cost matrix C. 

Proof. By definition of , the conditions 1. and 2. of Definition [23] are satisfied. From Theorem [2Tl it follows 
that properties 3. and 4. are also satisfied, so G^c a test family for MIPa,c- Minimality is due to the fact 
that removing any element {gfj,hfj) from G^c results in G^c \ {{9ij^ ^ij)}- However, this new set is not a test 
family since no oriented vector in G^c \ {idij^ ^ij)} '^^^ translated through a nonnegative vector in N" such 
that its tail meets a*. It is clear by definition that G^c depends only on A and C. □ 

Example 2.2. Let MIPa,c be the family of multiobjective problems, with the following constraints and objective 
function matrices: 

10 1 00" 

1 10 

Let (xi, X2, si, S2) be the vector of variables, where si and S2 are slack variables. In this example, using the 

order (see Remark\KB, Qc = {Gh, G^}, where Gl^^ = {"3 } = ((0, 1, 2, -1), (0, 1, 2, 0)), 

^1 = ((-1,1,0, -2),(0,1,0,0))} and ^ = = ((1, 0, 2, 0), (1, 0, 2, 0)), = ((1, -1, 0, 2), (1, 0, 0, 2))}. 

Figure\^ shows, on the [xi, X2) -plane, the -<c-skeleton of the fiber corresponding to the right-hand side vector 
(17, 11)*. In the box over the graph of the -<c-skeleton, we show the second components of the elements ofG^c- 
The reader may note that in the graph, the arrows have opposite directions due to the fact that the directed paths 
(improving solutions) are built subtracting the elements in G^c- describe how to compute the sets G\^. and 
G%^ in SectionWi 

Given G-ic i there are several ways to build a path from each feasible point in a fixed fiber to any Pareto-optimal 
solution. However, there is a canonical way to do it: Fix cr a permutation of the set {1, . . . , /i} and subtract 
from the initial point the elements of G"^^^^ for i = 1, . . . Add this element to an empty list. After each 
substraction by elements in GZ^^^) , i = l,...,Mi remove from the list those elements dominated by the new 
element. We prove in Section 3 that the result does not depend on the permutation a chosen. 

Example 2.2 (Continuation). This example shows the abovementioned different ways to compute paths from 
dominated solutions to any Pareto-optimal solution. The vector (9, 4, 9, 3) is a feasible solution for MIPa,c 
in the (17, iVf -fiber. Figure shows the sequence of Pareto-optimal points obtained from the feasible point 
(9,4,9,3) using the permutation ui = (1,2) (on the left) and using U2 = (2,1) (on the right). 



A 



2 2-10 
2 1 



C 
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9\ 



(4,5, 1 




(9, 0, 1, 11) 



Figure 2. The ^c-skeleton of the (17, ll)*-fiber of MIPa,c projected on the xi, a::2-plane. 




Figure 3. Two different ways to compute paths from (9, 4, 9, 3) to the Pareto-optimal solutions 
in its fiber. 

Remark 2.3. With the partial order, -<c, induced by C a directed path from a dominated point a to each 
Pareto-optimal point (3 in a fiber, applying the above method, cannot pass through any lattice point in this fiber 
more than fi times (recall that /x is the number of maximal chains in Q^^)- This implies that obtaining the 
Pareto-optimal solutions of a given MIPa^c using Q^^ cannot cycle. 

3. Test families and Partial Grobner bases 

In the previous section we motivate the importance of having a test family for MIPa,c since this structure allows 
us obtaining the entire set of Pareto-optimal solutions of the above family of multiobjective integer programs 
(when the right-hand side varies). Our goal in this section is to provide the necessary tools to construct test 
families for any multiobjective integer problem. Our construction builds upon an extension of Grobner bases 
on partial orders. 
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In order to introduee this structure we define the reduction of a pair {g, h) G Z" x Z" by a finite set of ordered 
pairs in Z" x 17^. Given is a coUection Qc Q Z" x Z!f: where Qc — {{gi,hi), . . . ,{gi,hi) : hk+i hk,k = 

The reduction of {g, h) by Qc consists of the process described in Algorithm [TJ 

Algorithm 1: Partial reduction algorithm 
input ■.R = {{g,h)},S = {{g,h)} 
For each {g, h) ^ S : repeat 

if h ~ gi and h — g are comparable by <c then 
I Ro = {{g - max^c{/i - g^, h - g})} 

else 

I Ro ^ {{g - gi,h- g^),{g - gi,h- g)} 
end 

For each r € Ro and s E R: 
if r -<c s then 
I R:^R\{sh 
end 

S Ro 
R :~ RU Ro', 
until {i : ft, - /i, > 0} = ; 

output: R, the partial reduction set of {g, h) by Qc 

The above reduction process extends to the case of a finite collection of ordered sets of pairs in Z" x Z" 
by establishing the sequence in which the sets of pairs arc considered. We denote by pRem{{g, h), {Qi))cr the 
reduction of the pair {g, h) by the family {Qi}l^i for a fixed sequence of indices a. The following result allows 
us to consider any sequence of indices for this process, since it establishes that the partial reduction does not 
depend on the chosen sequence. 

Theorem 3.1. Let Q be a finite set in Z" xZ" , whose maximal chains are Qi, ■ ■ ■ ,Qt, and a , a two permutations 
of the indices (1, . . . , t). Then, 

pRem{{g, h),Q)„ = pRem{{g, h),Q)y 

for each {g, h) eZ"" xZ^. 

Proof. Consider {g : g — g — ^ ^ Acr(i) sJ^'^Ij where Qi ~ {{gphi) : j = I . . . , ki}. It is clear that the 

i=i j=i 

elements in A^- does not depend on the permutation a, since reordering the sums does not give new elements. 

The elements in pRem{{g,h),Q)a are the element in Aq- deleting the comparable largest ones. Then, since 
^<7 = K' ' -P^em((5, h), Q)„ = pRem{{g, h), Q)^' . 

□ 

From now on. we denote by pRem{{g, h), Q) the set of remainders of {g, h) by the family Q = {5i}i=i 
sequence of indices. 

The reduction of a pair that represents a feasible solution, by a test family, gives the entire set of Pareto-optimal 
solutions. In order to obtain that test family, we introduce the notion of p-Grobner basis. This concept has 
been motivated by the fact that in the case where the ordering induced in N" by a single cost vector is total, 
a Grobner basis is a test set for the family of integer programs IPa,c (see [TU] or [iS] for extended details). In 
the single objective case the Buchberger algorithm computes the Grobner basis. However, in the multiobjective 
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case the cost matrix induces a partial order, so division or the Buchberger algorithm are not applicable. Using 
the above reduction algorithm we present in this section an "a la" Buchberger algorithm to compute the so 
called p-Grobner basis to solve MOILP problems. 

Definition 3.1 (Partial Grobncr basis). A family Q ~ {Gi, • ■ • , Gt} ^ Ia is a partial Grobner basis (p-Grobner 

basis) for the family of problems MIPa.c, if Git ■ ■ ,Gt are the maximal chains for the partially ordered set 
t 

(J Gi and for any {g,h) ^ Z" x Z'^; 

i=l 

g G Ker{A) pRem{{g, h),G) = {0}. 

A p-Grobner basis is said to be reduced if every element at each maximal chain cannot be obtained by reducing 
any other element of the same chain. 

Given a p-Grobner basis, computing a reduced p-Grobner basis is done by deleting the elements that can be 
reduced by other elements in the basis. After the removing process, the family is a p-Grobner basis having only 
non redundant elements. It is easy to see that the reduced p-Grobner basis for MIPa,c is unique and minimal, 
in the sense that no element can be removed from it maintaining the p-Grobner basis structure. 

This definition clearly extends to p-Grobner bases for the ideal Ia induced by A, once we fix the partial order, 
^(7, induced by C. 

The goal of this paper is to present algorithms to solve multiobjective problems analogous to the methods that 
solve the single objective case, using usual Grobner basis. These methods are based on computing the reduction 
of a feasible solution by the basis. The key for that result is the fact that the reduction of any pair of feasible 
solutions is the same, therefore the algorithm is valid for any initial feasible solution. The following lemma 
assures the same statement for the multiobjective case and p-Grobner bases. 

Lemma 3.1. Let G be the reduced p-Grobner basis for MIPa.c and ai,a2 two different feasible solutions in 
the same fiber of MIPa.c- Then, pRem{{ai, ai), G) = pRem{{a2, Oi2)^ G)- 

Proof. Let {(3,(3) £ pRem{{ai, ai) , G) , then (3 ~ a2 is in the same fiber and it cannot be reduced, so {(3,(3) S 
pRem{{a2,a2),G)- □ 

The following theorem states the relationship between the three structures introduced before: test families, 
reduced p-Grobner bases and the family G-<c- 

Theorem 3.2. The reduced p-Grobner basis for MIPa.c is the unique minimal test family for MIPa,c- 
Moreover, G~<cj introduced in is the reduced p-Grobner basis for MIPa^c- 

Proof. Let G — {Gi, • • ■ , Gt} be the reduced p-Grobner basis for MIPa,c- By definition of p-Grobner basis, it 
is clear that each Gi is totally ordered by its second component with respect to -<c (Condition 1). Condition 
2 follows because for each i and for each {g, h) £ Gi Q Z" x Z'l, clearly pRem{{g, h),G) — {0}, so 5 G Ker{A) 
and then A{h — g) = Ah. 

Now, let X G N" be a dominated solution for MIPA.c{b) then there is a Pareto-optimal solution, (3, such 
that (3 -<c X. By Lemma 13.11 pRem{{x, x),G) = pRem{{f3 , P) , G) , and by construction of the set of partial 
remainders, P G pRem{{P , P) , G) , and then x ^ pRem{{x,x),G)- This implies that {g,h) G Gi must exist such 
that X — gi -<c x, for some i G {1, . . . , t}. 

On the other hand, if x is a Pareto-optimal solution for MIPA,c{b), x G pRem{{x, x), G), and then, there exists 
no {g, h) in any Gi such that x ^ g -<c x. Therefore, for every i and for each {g, h) £ Gi, either x — g is infeasible 
or incomparable with x. 
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Minimality is due to the fact that removing an element from the reduced p-Grobner basis, that is the minimal 
partial Grobner basis that can be built for MIPa.c we cannot guarantee to have a test family because it may 
exist a pair (g, h) E Z" x Z" with g G Ker{A) that cannot be reduced to the zero set. 

The second statement of the theorem follows from Corollarv l2.2l □ 



In the following we describe an extended algorithm to compute a p-Grobner basis for Ia^ with respect to 
the partial order induced by C. First, for {g,h),{g , /i ) in Z" x Z" we denote by S^{{g,h),{g ,h )) and 
S^{{9,h.),{g' ,h')) the pairs 
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^9 


-2/i') 


if 7 - 


^9 


— 2h and 7 + 5 


— 2h are incomparable 



where 7 G N" and 7^ = max{/ii, /ij}, i ~ 1, . . . , n. 

The pairs S^{{g, h), {g ,h )) and S'^{{g, h), (g ,h )) are called 1 — Svector and 2 — Svector of {g, h) and {g ,h ), 
respectively. 

The notation is due to the analogy with the algebraic-geometrical notion of S-polynomial for a pair of polynomials 
with a given term order. Since we consider a partial order, it may happen that in the standard construction of a 
Svector (see |46j ) . we cannot decide which is the leading term. Therefore, in our definitions of Svectors we follow 
the standard construction but we must consider all possible combinations of leading terms, with respect to the 
partial order ^c- The following lemma is used in the proof of our extended criterion and it is an adaptation of 
the analogous result for total orders and usual S-polynomials. 

In the following, we denote by leadmond f) the set of leading monomials with respect to the order induced by 
C, for any multivariate polynomial / G K[xi,. . . ,a;„] and by 1-Spolynomial and 2-Spolynomial the binomial 
transcriptions of 1-Svector and 2-Svector (recall the equivalence between the pairs {u, v) and the binomial 
j,u-v — ii u is dominated by v). 



Lemma 3.2. Let /i, . . . , /s G K[xi, . . . , Xn] be such that there exists p G || leadmondfi) ■ Let f = Ci fi 

i=l i=l 

with Ci G K . If there exists q G leadmondf) such that q <c P, then f is a linear combination with coefficients 
in K of the k-Spolynomial, k = 1,2, of fi and fj, 1 < i < j < s. 



Proof. By hypothesis, /; ^ aip+ other smaller or incomparable terms, with a; G K, for all i. Then, / can be 

s s s 

rewritten as / = Ci fi = Ci Oip + other smaller or incomparable terms. Since q ~<c Pi then Ci Oi = 0. 

i—l i—1 i=l 
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By definition, for fc = 1, 2, = - ^ /j, tlius, 

/ = ci /i H h Cs /s 

= ci ai(— /i) H h Cs Z^) 

ai as 

= ci ai(— /i - — /2) + (ci fli + C2 a2) (— /2 - — /a) H 

fli a2 02 aa 

+ (ci ai H h Cs_i as-i) (^— /^-i -—/«) + (ci ai H 1- a^) — fs 

- rfi (/2,p)) + • • • + ds-1 S{{.fs-l,p), {fs,p)) + E ^ E ^((/-P)' 

i=l 'i=l 

i 

wfiere di = Cj aj . Tliis proves tlie lemma. □ 

i=i 

Tfie algorithm to compute standard Grobner bases is based on the Buchberger criterion, whose analogous for a 
partial order is the following. 

Theorem 3.3 (Extended Buchberger's criterion). Let Q — {^i, . . . ,Qt} with Qi C 1^ for alii = 1, . . . ,t, be the 

maximal chains of the partially ordered set {gi : gi G Gi, for some i — 1, . . . Then the following statements 
are equivalent: 

(1) Q is a p-Grdbner basis for the family MIPa.c- 

(2) For eachij = l,...,t and {g,h) G G,, {g',h') G Gj, pRem{S''iig,h),{g\h')),G) = {0} , for k =,1,2. 

Proof. The original Buchberger criterion was stated in a polynomial language. Therefore, we adapt our notation 
to follow the line of that proof. Each pair {u, v} is identified with the binomial — x^ , in the polynomial ring 
. . . ,Xn\, and our set Ia, with '^a = {x" ~ x^ : u — w G Ker{A)). The definition of partial remainders, 
pRem, is adapted accordingly. With these changes in the notation, the set setlm({u, v}) is identified with the 
elements in leadmonc{x^ — x"). 

Let 5 be a p-Grobner basis for i, j G {1, . . . , t} and {g,h) & Gi, {g ,h) G Gj- Then, S^{{g,K),{g ,h )), for 
A: = 1, 2, is in Ia, so by definition of p-Grobner basis, pRem{S^{{g, h), {g , h )), G) = {0}. 

Conversely, assume that for each (g, h) G Gi and {g ,h ) e Gj, pRem{S''{(g, h), {g , h )), G) 7^ {0} , for k =, 1, 2. 
Let (5, /i) G Z" X with g G Ker{A). We define f = x'^ - x^-^^ G Z[a;i, ...,Xn], and we denote by G* the 
polynomial set associated with G- 

Then, / can be written as a linear combinations of all the elements in G* (this representation is not unique): 

m 

i=l 

Let X = {Xi,...,X]^} be the set of maximal elements of the set {HiGi : Hi G leadmonc{hi),Gi G 
leadmoncigi)} , with respect to the partial order -<c'. 

If X = leadmoucif), the polynomial / can be partially reduced by the elements in G- This proves the result. 

Otherwise, assume that I G leadmonc{f)\X . Then, I comes from some simplification of the linear combination 
defining /. Then, the construction ensures that it must exist at least one element, Xi G X, such that I Xi. 

Set S = {i : HjGj = Xi with Hi G leadmonc{hi), Gi G leadmoncigi)}. For any j G S, we write hj = 
Hj + other terms and g = ^'^^Hjgj. Then, Xi G leadmonc{Hj gj), for all j G S. However, by hypothesis 

there exists G G leadmonc{g), with G ^ Xi. 
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Hence, by Lemma IX^ there exists d'l j. € K such that: 

2 

k=l r.s^S.r^s.Qs.gr&Gk 

Now, for any r, s G S*, = lc'm{Lr, L^) for some Lr G leadmonc{Hr gr) and G leadmondHs gs), so: 



15 



S'^{{Hr gr, Lr),{Hs gs, Ls)) = Y^iJrffr - H,, gs 

hr Lis 

where lr = h = and i/^.s lcm{lp{gr),lp{gs))- 

By hypothesis, pRem{S'^{gr, gs),G) = {0}. From the last equation we deduce that: 

pRemiS\Hrgr,Hsgs),g)^{0} 

this gives a representation: 

5'''(i^r 5r, 5s) = X! ^r,'^ 9y 

with g,y G G- 

max{iJ^^ G"': G leadmonc{h'^ s)t G'^ ^ leadmondgu)} = leadmonc{S{Hr gr, Hg gs)^) S, 



k 

r^s • 



By construction of S-polynomials, we have that there exists p G ^ such that p ^i, so, substitut- 
ing these expressions into g above and using that / = hj gj + hj gj = hj gj + g = hj gj + 

j^s jes j^s j^s 

dr.s <5'(i?s gs, Hr gr) = hj gj + h"^ g gy, we have expressed / as: 

i 

with one leading term, p, smaller than Xi. However, this is a contradiction proving the theorem. □ 



This criterion (the one in Theorem l3.3p allows us to describe a geometric algorithm which constructs a p-Grobner 
basis Qc for MIPa.c, and then a test family for that family of multiobjective problems. 

The first approach to compute a p-Grobner basis for a family of multiobjective programs, is an algorithm based 
on Conti and Traverso method for the single objective case [TU]. For this algorithm, the key is transforming the 
given multiobjective program into another one where computation is easier and an initial set of generators for 
I A are known. 

Notice that finding an initial set of generators for I a can be done by a straightforward modification of the Big-M 
method (see details, e.g. in [3]). 

Given the program MIPA.c{b), we consider the associate extended multiobjective program, EMIPA.c{b) as 

G C = {M ■ 1\C) G z("+l+")^^ Idm 





f 


-1 




the problem MIP^ ^ (b) where A = 






A 




\ 


-1 


J 



stands for the m x m identity matrix, AI is a large constant and 1 is the (m -I- 1) x fc matrix whose components 
are all 1. This problem adds m + 1 new variables, whose weights in the multiobjective function are big, and so, 
solving this extended minimization program allows us to solve directly the initial program MIPa.c- Indeed, 
any feasible solution to the original problem is a feasible solution to the extended problem with the first m 
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components equal to zero, so any feasible solution of the form (0, 0, ai, . . . ,a„) is non-dominated, upon 
the order -<g. by any solution without zeros in the first m components. Then, computing a p-Grobner basis for 
the extended program, allows us detecting infeasibility of the original problem. Furthermore, a trivial feasible 
solution, xo = (6i, . . . , 6m, 0, "t.^, 0), is known and the initial set of generators for Ia are given by {{Mi~Pi, Mi} : 
i = . . . , 7i} where Mi = (an — min{0, minj{aji}}, . . . , ami — min{0, imnj{aji}}, — min{0, rninj {aji}} , 0, 0), 
Pi = (0, 0|ei), for aU i = 1, . . . , 71, Mo = 1, 0, 0) and Pq = 0, Mi, P,, Mq, Pq € (see [2] 

for further details). 

Algorithm 2: Partial Buchberger algorithm I 

input : Fi = {Mo, A/i, . . . , M„} and F2 - {Pq, A, • ■ ■ , -Pn}, 

Mi ~ {an — min{0, minjlaji}}, . . . , ami — min{0, miiij{aji}}, — min{0, m.m.j{aji}}, 0, 0) (i > 0) 

Pi = {0, '"}+}, 0\ei) e (i > 0) 

A/o = (1, '"+.1,1,0,. ".,0) 

Po = (0,"+™.+i,0). 

repeat 

Compute, Qi, ■ ■ ■ ,Gt, the maximal chains for Q = (j){!F{Fi, F2)). 
for i, j € {1, . . . ,t}, i ^ j , and each pair {g, h) £ Gi, (g' , h') G Qj do 
Compute ^ pRem{S^{{g,h),{g' ,h!)),g), k = \,2. 

if R'' = {0} then 

I Continue with other pair. 

else 

I Add (j){T{r)) to Q, for each r £ R'' . 
end 
end 

until R^ = {0} for every pairs ; 

output: Q = {f/i, . . ■,Qq} 

p-Grobner basis for I a with respect to ^c- 

Then, we can state the following result. 

Theorem 3.4. Let G ^ {GiW^^ be a p-Grobner basis for EMI Pa^c- If {0,"}f},0, ai, . . . , a^) e pRem{{0,"}+} 
, 0, 61, ... , bn), G), then a = {ai, . . . , an) is a Pareto- optimal solution for MIPA.c{b). The entire set of Pareto- 
optimal solutions of MIPAfi{b) can be computed using the above construction. Moreover, if there are no a in 
the set pRem((0,b),G) whose m + 1 first components are zero MIPA,c{b) is infeasible. 

Proof. Let a be a vector obtained by successive reductions over G. It is clear that a is feasible because (0, a) is 
in the set of remainders of (0, /?) and then, in the same fiber. Besides, a is a Parcto-optimal solution because 
G is a. test family for the problem (Theorem 13. 2p . 

Now, if /?* is a Parcto-optimal solution, by Lemma [3.11 pRem (( 3* . 3*).G)) ~ pRem{{P , /3) , G)) , but since (3* is 
a Parcto-optimal solution, it cannot be reduced so (/?*,/?*) G pRem{{(3* , i3*),G)), and then, also to the list of 
partial remainders of {(3,f3) by G- O 

Hosten and Sturmfels [26j improved the method by Conti and Traverso to solve single-objective programs 
using standard Grobner bases. Their improvement comes from the fact that it is not necessary to increase 
the number of variables in the problem, as Conti and Traverso's algorithm does. Hosten-Sturmfels's algorithm 
allows decreasing the number of steps in the computation of the Grobner basis, but on the other hand, it needs 
an algorithm to compute an initial feasible solution, that in Conti and Traverso algorithm was trivial. We have 
modified this alternative algorithm to be used to compute the entire set of Pareto-optimal solutions. The first 
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step in the algorithm is computing an initial basis for the polynomial toric ideal S^i = {x^ —x"" : u — v G Ker{A)), 
that we are identifying with I a. This step does not depend on the order induced by the objective function, 
so it can be used to solve multiobjectivc problems. Details can be seen in |26j . Algorithm [3] implements 
the computation of the set of generators of 3a- This procedure uses the notion of LLL-reduced basis (see 
|32| for further details). In addition, we use a w-graded reverse lexicographic term order, ~<^^\ induced by 
Xi+i > ■ ■ ■ > Xi-i > Xi (with Xn+i := xi), that is defined as follows: 



\ Y.7=i '^iO^t = Y.t=i ^iPi and a <iex P 
where lo G W] is chosen such that Xi+i > • • • > Xi^i > xi. 



Algorithm 3; setof generators (A) 



input : A e Z™^" 

(1) Find a lattice basis B for Ker{A) (using the Hermite Normal Form). 

(2) Replace B by the LLL-reduced lattice basis Bred in the sense of Lovasz (see [35] for more details). 

Let Jo (a;"+ - x"- : li G Bred)- 
for i — 1, . . . , n do 

Compute Ji ~ (Ji_i : x°°) as: 

(a) Compute Qi-i the reduced Grobner basis for Ji-i with respect to ■ 

(b) Divide each element / G Qi-i by the highest power of Xi that divides /. 
end 

output: 3^ := J„ = {x^^ — .t"i , . . . , — x""" } system of generators for Ia- 



'^A consists of binomials x"' — x^' with Ui — Vi € Ker(A), for i = 1, . . . , s. Coming back to our notation, each 
binomial, x" — x^, in 3^ is identified with {u, v} € Ia, so computing a set of generators for 3^ gives us, in some 
sense, a finite number of generators for the set that represents the constraints matrix. We compute in the next 
step a partial Grobner basis from the initial sets Fi = {wi, . . . ,Us} and F2 = {vi, . . . ,Vs} using our extended 
Buchberger algorithm: 

Algorithm 4: pgrobner (Fi, F2) 

input : Fi {Mi, . . . , Mr} and F2 = {Pi, . . . , P,.}- 

repeat 

Compute, Gi, ■ ■ ■ ,Gt, the maximal chains for G = (f){J^{Fi, F2)). 
for i, j G {1, . . . , t}, i ^ j, and each pair (g, h) G Gi, (5', h') G Gj do 
Compute ^ pRem{S''{{g,h),{g' ,h')),G), k = 1,2. 

if R'' {0} then 

I Continue with other pair. 

else 

I Add ^(.F(r)) to G, for each r e R''. 
end 
end 

until R'' = {0} for every pairs ; 
output: G = {Gi, ■ ■ - jGq} 

p-Grobner basis for the set spanned by {{Mi, Pi} : i = 1, . . . ,r} with respect to ^c- 
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Once wG have obtained the partial Grobner basis using the above algorithm, we can compute the entire set of 
Pareto-optimal solutions for MIPA,c(b) by the following algorithm: 

Algorithm 5: Pareto-optimal solutions computation for MIPA,c{b) 
input : MIPA.cib) 

Step 1. : Compute an initial feasible solution, ag, for MIPA,c{b)- It consists of finding a solution for 

the diophantine system of equations Ax = 6, a; G Z". 
Step 2. : Compute a system of generators for I a'- {{ui, Vi} : i = 1, . . . , s}, using setof generators {A) . 
Step 3. : Compute the partial reduced Grobner basis for MIPa.c, 0c = {Qi, ■ ■ ■ , Gt}, using pgrobner ( 

Fi, F2 ) , where Fi = {ui : i ~ 1, . . . ,r} and F2 ~ {vi : i ~ 1, . . . ,r}. 
Step 4. : Calculate the set of partial remainders: R := pRem{ao,Gc)- 
output: Pareto-optimal Solutions : R. 

There are some interesting cases where our methodology is highly simplified due to the structure of the set 
of constraints. One of these cases is when the dimension of the set of constraints is n — 1. The next remark 
explains how the algorithm simplifies in this case. 

Remark 3.1. Let A he a m x n integer matrix with rank n — 1. Then, since dim{Ker{A)) ~ 1, the system 
of generators for I a (Step 2) has just one element, {g,h), and the p-Grdbner basis (Step 3) is the family 
Q = {{(.9,/i)}} because no Svector appears during the computation of the Buchberger algorithm. In this case, 
Pareto-optimal solutions are obtained as partial remainders of an initial feasible solution {a, a) by {g,h), i.e., 
the entire set of Pareto-optimal solutions is a subset of T = {a — Xg : X E Z+}. More explicitly, the set of 
Pareto-optimal solutions for MIPA.c{b) is the set of minimal elements (with respect to <c) ofT. 

In order to illustrate the above algorithm, we present an example of MOILP with two objectives where all the 
computations arc done in detail. 

Example 3.1. 

min {lOx -\- y , x -\- lOy} 
s.a. 

2x + 2y ^17 
2y 11 

x 10 

x,y e z+ 

Transforming the problem to the standard form results in: 



min {10.T + y + 0z + 0t + Qq,x+ lOy + Oz + Qt + Oq} 
s.a. 

2x-\-2y- z =17 
2y + t =11 
x + q =10 
x,y,z,t,q e Z+ 

Step 1. : Feasible solution for MI PA,c{b): w= (9,4,9,3,1). 
Step 2. : Following the steps of Algorithm\3i 

(1) Basis for Ker {A) ; B := {(0, 1, 2, -2, 0), (-1, 0, -2, 0, 1)}. 

(2) LLL-reduced basis for B : Bred B := {(-1, 0, -2, 0, 1), (-1, 1, 0, -2, 1)}. 

(3) Jo := (a;"+ - a;"- : u e Bred) = (2^5 - xixl,X2X^ - xixD 
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Ia - ({((1,0,0, 0,1), (0,1, 0,2,0)), ((1,0, 2,0,0), (0,0,0, 0,1)), ((0,1, 2, 0,0), (0,0, 0,2,0))} 
Step 3. : Computing a p-Grdbner basis for I A, using the order ^'^(^ (Remark \2.2\) . and following Algorithm 
[7] we obtain G , whose maximal chains are: 

Gi: {((0, 1, 2, 0, 0), (0, 0, 0, 2, 0), (0, 1, 2, 0, 0)), ((0, 1, 0, 0, 2), (2, 0, 2, 2, 0), (0, 1, 0, 0, 2)), 

((0, 1,0, 0,1),(1,0, 0,2,0), (0,1,0,0,1))}. 
G2: {((1,0, 0,4,0), (0,2,2, 0,1), (1,0, 0,4,0)), ((1,0,2,0,0), (0,0,0, 0,1), (1,0,2,0,0)), 
((1,0,0,2,0),(0,1,0,0,1), (1,0,0,2,0))}. 
Step 4. : Partial remainders: Reducing first by Gi •' 
pi?em((9,4,9,3,l),ai) - {(9,0,1,11,1)}. 
Then, reducing each remainder by G2-' 

pi?em((9, 0, 1, 11, 1), 02) = {(9, 0, 1, 11, 1), (8, 2, 3, 7, 2), (7, 2, 1, 9, 3), (6, 3, 1, 5, 4), (5, 4, 1, 3, 5), (4, 5, 1, 1, 6)}. 

The entire set of Pareto- optimal solutions is: 

{(9,0, 1,11,1), (8, 1,1, 9, 2), (7, 2, 1,7, 3), (6,3, 1,5,4), (5,4, 1,3, 5), (4,5, 1,1, 6)} 

Figure^ shows the feasible region and the Pareto- optimal solutions of the example above. 

y 6' 



5 



4 

















































1 c 

. . . ., 


') c 


') < 


') < 


3 < 


3 <) 












1 c 

■V •« 


3 c 


') < 


5 C 


3 ( ) 














1 { 

■V < 


3 c 


5 < 


> <) 






I 










1 c 


D C 

» C 


3 ( ) 




















> ( ) 

►t © 



0123456789 10 

X 



Figure 4. Feasible region, Pareto-optimal solutions and improvement cone for Example 13. II 



In addition, we have evaluated the problem with the same feasible region but choosing a cost matrix such that the 
respective normal vectors of each of the rows in the matrix form an acute angle. Then, non supported solutions 
appear in the set of Pareto-optimal solutions. Figure\Bi shows the Pareto-optimal solutions for the same feasible 
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Figure 5. Feasible region, Parcto-optimal solutions and improvement cone for Example 13.11 
withC= [[10, -!],[-!, 10]] 

4. Computational Results 

A series of computational experiments have been performed in order to evaluate the behavior of the proposed 
solution method. Programs have been coded in MAPLE 10 and executed in a PC with an Intel Pentium 
4 processor at 2.66GHz and 1 GB of RAM. In the implementation of Algorithm |4] to obtain the p-Grobner 
basis, the package poset for Maple [IHI has been used to compute, at each iteration, the maximal chains for 
the p-Grobner basis. The implementation has been done in a symbolic programming language, available upon 
request, in order to make the access easy to both optimizers and algebraic geometers. 

The performance of the algorithm was tested on randomly generated instances for knapsack and transportation 
multiobjective problems for 2, 3 and 4 objectives. For the knapsack problems, 4, 5 and 6 variables programs 
were considered, and for each group, the coefhcients of the constraint were randomly generated in [0, 20] and 
the coefficients of the objective matrices range in [0, 20]. Once the constraint vector, (ai, . . . , a„), is generated, 
the right-hand side is fixed as 6 = J27=i ensure feasibility. 

The computational tests have been done in the following way for each number of variables: (1) Generate 5 
constraint vectors and compute the initial system of generators for each of them using Algorithm^ (2) Generate 
five random objective matrices for each number of objectives (2, 3 and 4) and compute the corresponding p- 
Grobner basis using Algorithm |31 and (3) with h ~ \^ X]"=i \ ^'^^ ^^r each objective matrix, compute the 
Pareto-optimal solutions using Algorithm [5l 

Table[l]contains a summary of the average results obtained for the considered knapsack multiobjective problems. 
The second, third and fourth columns show the average CPU times for each stage in the algorithm: sog is the 
CPU time for computing the system of generators, pgrobner is the CPU time for computing a p-Grobner basis, 
and pos is the time for computing a feasible solution and partially reduce it to obtain the set of Pareto-optimal 
solutions. The fifth column shows the total time for computing the set of Pareto-optimal solutions for the 
problem. Finally, the sixth and seventh columns show the average number of Pareto-optimal solutions and the 
number of maximal chains in the p-Grobner basis for the problem. The problems have been named as knapN_0 
where N is the number of variables and is the number of objectives. 

For the transportation problems, instances with 3 origins x 2 destinations, 3 origins x 3 destinations and 4 
origins x 2 destinations were considered. In this case, for fixed numbers of origins, s, and destinations, d, the 
constraint matrix, A € ^(s+'^)x(«'')^ is fixed. Then, we have generated 5 instances for each problem of size s x d. 
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Table 1. Summary of computational experiments for knapsack problems 



Each of these instances is combined with 5 different right-hand side vectors. The procedure is analogous to 
the knapsack computational test: a first step where a system of generators is computed, a second one, where 
the p-Grobner basis is built and in the last step, the set of Parcto-optimal solutions is computed using partial 
reductions. Table [H shows the average CPU times and the average number of Pareto-optimal solutions and 
maximal chains in the p-Grobner basis for each problem. The step column shows the average number of steps 
in the p-Grobner computation, and act_pGB is the average CPU time in the computation of the p-Grobncr 
basis elapsed since the last element was added to the basis until the end of the process. The problems have been 
named as transNxM_0 where N is the number of origins, M is the number of destinations and is the number of 
objectives. 
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Table 2. Summary of computational experiments for the battery of multiobjective transporta- 
tion problems 



As can be seen in tables [T] and O the overall CPU times are clearly divided into the three steps, being the 
most costly the computation of the p-Grobner basis. In all the cases more than 99% of the total time is spent 
computing the p-Grobner basis. Once this structure is computed, obtaining the Pareto-optimal solutions is 
done very efficiently. 

The CPU times and sizes in the different steps of the algorithm arc highly sensitive to the number of variables. 
However, our algorithm is not very sensitive to the number of objectives, since the increment of CPU times 
with respect to the number of objectives is much smaller than the one with respect to the number of variables. 

It is clear that one can not expect fast algorithms for solving MOILP, since all these problems are NP-hard. 
Nevertheless, our approach gives exact tools that apart from solving these problems, give insights into the 
geometric and algebraic nature of the problem. 

As mention above, using our methodology one can identify the common algebraic structure within any multi- 
objective integer linear problem. This connection allows to improve the efficiency of our algorithm making use 
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of any advance that improves the computation of Grobner bases. In fact, any improvements of the standard 
Grobner bases theory may have an impact in improving the performance of this algorithm. In particular, one 
can expect improvements in the efficiency of our algorithm based on the special structure of the integer program 
(see for instance Remark l3.ip . In addition, we have to mention another important issue in our methodology. As 
shown in Theorem l3.21 solving MOILP with the same constraint and objective matrices requires computing only 
once the p-Grobner basis. Therefore, once this is done, we can solve different instances varying the right-hand 
side very quickly. 

Finally, we have observed from our computational tests that a significant amount of the time, more than 60% 
of the time (see column act_pGB), for the computation of the p-Grobner basis is spent checking that no new 
elements are needed in this structure. This implies that the actual p-Grobner basis is obtained much earlier 
than when the final test is finished. A different truncation strategy may be based on the number of steps 
required to obtain the p-Grobner basis. According to the exact method, the algorithm stops once in a step no 
new elements are added to the structure. Our tables show that in most cases the number of steps is 2, actually 
only one step is required to generate the entire p-Grobner basis (see column steps). These facts can be used to 
accelerate the computational times at the price of obtaining only heuristic Pareto-optimal solutions. This idea 
may be considered an alternative primal heuristic in MOILP and will be the subject of further research. 
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